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Abstract 



We analyze the capability of prompt photon production in pp and pp collisions to con- 
strain the gluon distribution of the proton, considering data from fixed-target experiments 
as well as collider measurements. Combined fits are performed to these large-pr direct-7 
cross sections and lepton-proton deep-inelastic scattering data in the framework of next- 
to-leading order perturbative QCD. Special attention is paid to theoretical uncertainties 
originating from the scale dependence of the results and from the fragmentation contri- 
bution to the prompt photon cross section. 



* On leave of absence from Sektion Physik, Universitat Munchen, D-80333 Munich, Germany 



1 Introduction 



The production of high-p^ prompt photons in pp or pp collisions provides an important 
probe of the proton's gluon distribution, g(x, Q 2 ), due to the presence and dominance of 
the leading order (LO) 0(aa s ) 'Compton-like' subprocess qg — > jq. In fact, constraints 
on g(x, Q 2 ) for 0.3 ~ x ~ 0.6 at Q 2 ~ 10 GeV 2 , derived mainly from the WA70 fixed-target 
pp — > data have been the backbone of the gluon determination in many parton 
density analyses J|, [3], f|, [5| ever since the pioneering work of ||. In that paper, a combined 
next-to-leading order (NLO) fit to the WA70 data and to deep-inelastic scattering (DIS) 
results from the BCDMS collaboration |7j was performed, and available direct photon data 
from other fixed-target experiments || |D| as well as from ISR |T0L |TTJ and SppS [|T2|, |T3] ] 
were compared to this fit. Since then, major theoretical and experimental developments 
concerning direct-7 production have taken place, and HERA results begin to add DIS 
constraints in a previously unexplored kinematic region. Therefore, we feel that it is 
time now for a reanalysis of the prompt photon data and their implications on the gluon 
density. 



Experimentally, much progress has been made during the last years. Some of the direct 
photon data sets mentioned above have been superseded by improved analyses [14], [1^, [16|] , 
usually providing smaller statistical and systematic errors. Even more importantly, the 
accessible range of fractional gluon momenta x has been considerably enlarged by the 
recent, partly very precise measurements at the Fermilab Tevatron |17|, [Tj| pjj. By now, 



the region 0.01 ~ x ~ 0.6 is completely covered, with the data from WA70 [TJ], R806 
n I , UA2 [[HJ and CDF presently dominating in their respective kinematical regimes. 



Also, the DIS kinematic coverage has been dramatically improved by the NMC and 
recent HERA |21], [EJ F^{x, Q 2 ) structure function data. Quark density measurements, 



extending down to almost x = 10~ 4 now, also imply an important constraint on the gluon 
distribution due to the momentum sum rule. Moreover, scaling violations of F 2 P at HERA 
now begin to constrain g(x, Q 2 ) severely, especially at x ~ 0.01 |23[| . Hence it is interesting 



to examine quantitatively the questions of whether a successful NLO perturbative QCD 
description of all these data is possible and how much freedom is left for g(x, Q 2 ). 

On the theoretical side, it was not yet possible to perform a complete and fully con- 
sistent calculation of the NLO prompt photon cross section at the time of the analysis in 
0, since the NLO fragmentation contribution, based on the partonic 2^3 QCD subpro- 
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cesses and on corresponding parton-to-photon fragmentation functions |25|, [26|| , was 
not yet available. Also, the development of a proper NLO theoretical implementation of 
isolation cuts imposed on the cross section in the high-energy SppS and Tevatron 

experiments, was only recently finished [ 29] and demonstrated to be phenomenologically 
important |3T| . Even in a very recent global study on direct-7 production by members of 



the CTEQ group |3l| , the fragmentation pieces are included only in LO and the isolated 
collider data are transformed to fully inclusive cross sections on this basis. 



Finally, the theoretical uncertainties still present in the NLO treatment have not been 
considered to their full extent in the literature so far. For instance, the dependence of 
the theoretical cross section on unphysical scales, such as the renormalization scale [ir 
and the factorization scale /if, has been treated in the fits of || |32| by optimizing the 
scales using the 'principle of minimal sensitivity' EM, i.e., by choosing scales where the 



NLO cross section is stationary with respect to small changes in /i^ and /ip. While this 
concept certainly is attractive to some extent, see however [J27I] for a critical discussion, 
it inevitably does suppress the uncertainty of the theoretical prediction due to the scale 
dependences. The recent study of |3l|] has not thoroughly addressed this issue either: in 
that analysis the scales were included, subject to the constraint /i# = /i^, among the 
parameters fitted to data, which represents just another kind of optimization. There are 
also other ambiguities, e.g. originating from the experimentally virtually unknown parton- 
to-photon fragmentation functions. In order to examine the question of how effective the 
constraints on the gluon distribution coming from prompt photon production really are, 
it is a crucial issue to take into account such uncertainties inherent to the calculation. 



The remainder of this paper is organized as follows: In section 2 we briefly recall 
the main ingredients needed to calculate the inclusive and isolated prompt photon cross 
sections in NLO. Section 3 is devoted to a detailed discussion of the main theoretical 
uncertainties mentioned above and their effects on the calculated cross sections. In section 
4 we present combined analyses of an exhaustive set of DIS and prompt photon data in 
order to arrive at conclusions about the quantitative viability of the NLO framework and 
about the constraints on the gluon distribution. Finally, we summarize our main findings 
in section 5. 
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2 General Framework 



Two types of processes contribute to the prompt photon production cross section: the 
so-called 'direct' piece, where the photon is emitted via a pointlike (direct) coupling to a 
quark, and the fragmentation piece, in which the photon originates from the fragmentation 
of a final state parton. Despite the fact that its corresponding partonic subprocesses are 
of order a 2 , the fragmentation contribution is present already in LO since the parton-to- 
photon fragmentation functions are effectively of order a/a s in perturbative QCD, where a 
denotes the fine structure constant. Next order corrections in the strong coupling constant 
a s have been calculated in the MS renormalization and factorization schemes for both the 
direct [T2^, [J3], and the fragmentation subprocesses, hence the cross sections can be 
consistently calculated to order aa 2 s . The cross section for the fully inclusive production 
of a prompt photon with momentum p 7 schematically reads 

da = da dir + da frag = J2 J dx a dx b f a (x a , fi 2 F )f b (x b , fx F ) x (1) 

a,b=q,q,g 

—da c ab (p-y, x a , x b , z, n Rl h f , M f )D2(z, M 2 f ) 

where z m i n = xt cosh rj with the prompt photon's rapidity rj, and xt = 2pp/^. In 
eq. (1), da l ab represent the subprocess cross sections for partons a, b producing a particle 
i {i — 7, q, g), integrated over the full phase space of all other final state particles. 
fi(x, fip) denotes the number density of the parton type i in the proton (or antiproton) at 
momentum fraction x and scale /i^, and D^(z, Mp) is the photon fragmentation function 
at scale Mp, z being the fraction of energy of the fragmenting parton c transferred to the 
photon. 

As already mentioned in the introduction and made explicit in eq. (1), the cross section 
in any fixed order of perturbation theory depends on unphysical scales which have to be 
introduced in the procedure of renormalization (fin) and of factorization of initial (i^f) 
and final (Mp) state mass singularities. The latter type of singularities appears, e.g., in 
the calculation of the 0(a^) ab — > cde NLO fragmentation subprocess cross sections when 
the parton c (which then fragments into the photon) becomes collinear with particle 
d or e, but also in the calculation of the 0(aa 2 ) ab — > 7<ie NLO 'direct' subprocess 
cross sections when the photon and a final state quark are collinear. These singularities 
need to be factorized at a scale Mp into the 'bare' fragmentation functions in order to 
render the cross section finite. The fragmentation functions then obey corresponding 
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NLO evolution equations. Since factorizing singularities is not a unique procedure but 
depends on the factorization prescription adopted, it becomes obvious that only the sum 
of the direct and the fragmentation pieces is a physical (scheme independent) quantity 
beyond the LO, but not these parts individually. In particular, even if the fragmentation 
contribution turns out to be numerically small, its addition on a LO basis to a NLO 'direct' 
piece is theoretically inconsistent and yields a scheme dependent cross section. Needless 
to say that a consistent NLO calculation also affords parton distributions and photon 
fragmentation functions evolved according to their respective NLO evolution equations. 
For calculating these (Revolutions we use the Mellin-n space technique described in 



At very high-energy pp colliders the photon is experimentally required to be 'isolated' 
in order to suppress the huge background due to copious production of n°, the decay of 
which can fake a prompt photon event. This means that the amount of hadronic energy 



E had allowed in a cone y(A0) 2 + (A77) 2 < R around the photon direction is limited to 
a small fraction of the photon energy, E^ad < tE~ with e ~ 0.1. In order to compare 
QCD predictions with isolated collider data, the theoretical calculation has to include 
this isolation criterion which leads to a significant decrease from the fully inclusive cross 



section |29|, |37|. In |29fl a simple, yet accurate way of incorporating the isolation cut into 
the NLO calculation has been developed. Starting from the fully inclusive cross section, 
the isolated one is obtained by introducing a subtraction term |27|, p9[ , 

da isol {R, e) = da md - da sub (R, e) , (2) 

where da sub (R, e) is the cross section for having more hadronic energy than eE^ in the 
cone around the photon. It turns out that da sub (R, e) can be easily and reliably calculated 
in the approximation of a rather narrow cone. An equation similar to eq. (2) (which holds 
for the NLO direct contributions) can also be written down for the NLO fragmentation 



piece |29| , in which case the main effect of imposing the isolation cut is to raise the lower 
integration limit in eq. (1) to z min = 1/(1 + e). In this way, it is possible to calculate also 
the isolated prompt photon production cross section measured at high-energy colliders in 
a consistent way beyond the leading order. 

We conclude this section by noting that the recent study in |U| reports a (rather 



Pt independent) 0(10%) discrepancy between the NLO programs of p8| and |35|] in the 



kinematic regime of the CDF measurements, which persists also if the fully inclusive cross 



section is considered |38[ . In this context it is interesting to mention that there is no such 



discrepancy between the two calculations of [33] and [35|, which are in exact agreement 
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concerning the NLO direct contribution to the fully inclusive cross section for all pr and 



■y/s. For all calculations to follow we use the program of [35| for the direct part of the 
NLO prompt photon cross section, along with the expressions of pl|| for the fragmentation 
contribution. When calculating the isolated cross section we complement these programs 
according to the prescription of 



3 Theoretical Uncertainties 



In this section, we address the main uncertainties entering the NLO calculation of the 
prompt photon cross section, namely the dependence on the photon fragmentation func- 
tions and on the renormalization and mass factorization scales. As a point of refer- 
ence, we first calculate the cross section for a fixed 'standard' set of input distributions 
and parameters and confront it with the data. For this purpose, the parton distribu- 



tions and photon fragmentation functions are taken from GRV j|, 26fl , together with the 



value of the QCD scale parameter for four active flavours, = 200 MeV. We choose 
fJ>R — f^F — Mp = pt/2 for the renormalization and factorization scales except for the 
isolated prompt photon data, where Mp = Rpp seems more appropriate |2"7fl . 



The very small charm effects in the cross section at fixed-target and ISR energies 
are neglected. For the SppS and the Tevatron experiments, however, charm-induced 
contributions coming, e.g., from eg —>■ are not negligible. We employ the effective 
(massless) charm quark distribution of |§ in calculating these contributions. We have 
checked that reasonable variations of this charm density do not significantly alter our 
results. An alternative approach to the heavy quark (charm) contribution is to perform 
mass factorization only for the light u, d, and s quarks. In this scheme, also used in 
0], the heavy flavours c, b, ... do not act as partons in the proton, and in LO charmed 
prompt photon events are only introduced via the processes gg — *> 'jee and qq — > ^cc with 



massive charm quarks. The results of |39| show, however, that these two approaches yield 
very similar results at least in LO. Thus the theoretical uncertainty originating from the 
charm treatment seems to be rather small. For the rest of this paper, we therefore use 
the fixed intrinsic charm quark distribution of || which facilitates the calculations. 

The sets of experimental data we take into account in this section are the fixed-target 
data of WA70 Q (y/s = 24 GeV, pp inclusive), the ISR results from R806/7 @ (y/s = 63 
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GeV, pp inclusive), the SppS results from UA2 |TJj (y^s = 630 GeV, pp isolated) and the 



Tevatron data of CDF @ (y/H = 1.8 TeV, pp isolated with R = 0.7 and e = 2 GeV/p T )- 
As mentioned above, these data sets dominate in their respective kinematical domains. 



The cross sections of W and [18 1 have been averaged over the experimentally covered 



regions of rapidity 77, while the results in \TA\ and [T5|] have been presented at 77 = 0. 
For our comparisons and fits, we add the statistical and systematic errors in quadrature, 
using point-to-point errors where these are separately available |15|, [[HQ- Fig. 1 displays 



the results for our 'standard' choice of input distributions and parameters. We show the 
'default quantity' (data — theory) /theory versus xt- This provides a particularly easy 
visualization of the (dis) agreement between data and theoretical calculation in view of 
the strong px fall-off of the cross section, xt is a good representative of the Bj0rken-x 
values predominantly probed in the gluon distribution at given pt and y/s. 

As can be seen from Fig. 1, the overall agreement between data and the NLO theoret- 
ical prediction is good though not complete, especially if the very small errors of the CDF 
data are taken at face value. The agreement between theory and the fixed-target and 
ISR measurements is very good, whereas the comparison with the high-energy collider 
results seems to be slightly less successful, since both the CDF and the UA2 data show 
a somewhat stronger rise for small xt than the theoretical cross sectionP]. This effect is 
more pronounced and statistically more meaningful for the CDF results, which possess 
the smallest point-to-point errors of all data sets and therefore provide a very precise 
measurement of the slope of the cross section. Note that the CDF [T3 as well as the 



UA2 [15] data are subject to a normalization uncertainty of about 10% which we have 
used in the figure to center the results on the zero line. It has to be emphasized that a 
much stronger discrepancy between high-energy collider data and NLO calculations was 
reported previously ||T7| , [18| . As was shown in | 30| , a dramatic improvement is obtained 



by using 'modern' sets of (steep) parton distribution functions like, e.g., those of GRV 
||. [|] or the most recent MRS(A',G) sets [§] as well as, equally important, by includ- 
ing a properly isolated NLO fragmentation contribution in the calculation. According to 
Fig. 1, the latter amounts to a 20% slope effect for CDF conditions, thus its inclusion 
is clearly crucial for a quantitative comparison between the experimental and theoretical 
cross sections. The fragmentation piece is non-negligible despite the presence of the iso- 

1 We thank S. Kuhlman for providing the break-up of the systematic errors of the CDF direct photon 
data into a point-to-point and a fully correlated part. 

2 Note, however, that there is good agreement between NLO theory and the preliminary Tevatron DO 
results |19) , which on the other hand have sizeably larger errors than the CDF data . 



6 



lation cut without which it would easily contribute about 50% to the total cross section 
f25| , |30| , j37| . The results in Fig. 1 indicate a (minor) remaining discrepancy between the 
PT-slopes of the experimental and theoretical cross sections even after the improvements 
of |30| have been applied. One furthermore infers from the figure that fragmentation also 



plays an important role in the calculations in the ISR and fixed-target regions, where no 
isolation cut has been applied. Here it leads to an effect of partly even more than 20%, 
but influences the slope to a lesser extent. 

An important uncertainty in the calculation is the dependence of the cross section 
on the parton-to-photon fragmentation functions which are experimentally unknown so 
far. Two partly very different NLO sets of such distributions have been suggested in the 
literature, namely in |25j (ACFGP) and in [|(| (GRV). Fig. 2 presents and DJ from 
both groups at two scales relevant for our comparisons to data. In both [25[] and [26[ . 



the fragmentation functions are assumed to evolve from a pure vector meson dominance 
(VMD) input at some very low scale. However, this boundary condition for D^(z,Q 2 ) 
has been implemented in rather different factorization schemes. ACFGP use the MS 
scheme, whereas GRV impose the VMD input in the timelike version of the so-called 
DIS 7 scheme, originally introduced for the (spacelike) parton structure of the photon f40| . 
In the MS scheme employed here the latter ansatz corresponds to an additional, rather 



large input for D^, which guarantees the positivity of the timelike structure function /j 



(T) 



for single photon inclusive e + e annihilation, e + e — > jX. Thus the quark-to-photon 
fragmentation functions of GRV [EBJ are larger than the ones of ACFGP [|25[|, especially 



at low scales, despite the fact that in |25[ a sizeably larger VMD input is employed. 
On the other hand, DJ of ACFGP is much larger than its GRV counterpart. This 

) crlnnnir innnt in hasprl on thp ascnimntinn T,^ _ 



is due to the huge VMD gluonic input in p5[, based on the assumption Df = D£° 



with Jq dz zDg° (z,2 GeV 2 ) = 0.5, meaning that as much as half of an outgoing gluon's 
momentum is carried away by neutral pions alone. 

The effects of these differences are displayed in Fig. 3, where we show the (scale and 
scheme dependent) relative importance of the fragmentation part of the prompt photon 
cross sections for ISR and Tevatron energies both for the total and for the gluon initiated 
contributions. We have normalized all results to the direct cross section, and apart from 
the fragmentation functions, all parameters and distributions have been chosen as in 



the 'standard' calculation above. As already obvious from Fig. 1, the ACFGP |25j and 
GRV |26[ fragmentation functions yield very similar results for the total fragmentation 
contribution. For the isolated Tevatron case, g — > 7 fragmentation plays an almost 
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negligible role due to the high z m i n cut implied by the isolation criterion, e.g, z m i n — (1 + 
2 GeV/pr) -1 for CDF conditions. The q —>■ 7 pieces are rather similar, since D^(z,Q 2 ) 
is probed at large z and Q 2 here. In the fully inclusive ISR case, on the other hand, 



z min = xt = 2pt/\/s ~ 0.15 for 7] = 0, and the huge difference in DJ between |25j and 
|213fl enters the fragmentation cross section. However, since the scales Q 2 ~ Vr/^ are rather 
low here, this effect is strongly compensated by the difference in the quark contributions. 
Of course, the difference in the total results is not necessarily fully representative of the 
uncertainty originating in the fragmentation part of the cross section. We have checked 
that in the theoretically more realistic GRV |26| 50% change in the VMD input 



distributions for D? has no sizeable effect on the results due to the dominance of D' 



Hence we will keep the fragmentation functions of |26| for the rest of this paper. Clearly, 
experimental information on the fragmentation functions, e.g. from e + e~ — > 7X jHJ, is 
needed. 

Let us now discuss the scale dependence of the results, i.e. the changes in the theoretical 
predictions for varying fi R and fip. It turns out that the NLO cross section depends only 
very weakly on Mp both for the isolated |29| and the fully inclusive cases. For the latter, 



e.g., the difference of the results for Mp = pp/2 and Mp = pp at fixed //r, /if does 
not exceed 1%. We therefore keep the fragmentation scale Mp fixed at Mp = pp/2 for 
the fixed-target and ISR experiments and Mp = Rp T for the isolated cross sections in 
the following. The dashed and the dash-dotted lines in Fig. 4 display the shifts in the 
theoretical results if we choose /i# = \ip = 0.3 pp or [ir = fip = 1.0 pp, respectively, 
instead of [Ir = \ip = pp/2. More precisely, the curves show (cr' th — cr th )/a th , where a' th 
is the theoretical cross section as calculated with the new values for the scales, whereas 
Oth corresponds to the 'standard' calculation. It becomes obvious that the results for 
fJ'R = = 0.3 pp or 1.0 pp amount to almost a constant shift in the normalization of 
the theoretical cross section as far as the CDF and UA2 data are concerned, and do not 
provide a change in the slope of the cross section. It can also be seen that the theoretical 
cross section at lower energies shows a rather strong scale dependence. So far, our results 



are in agreement with the claims in [|18|, [31] that scale uncertainties provide (almost) no 



effect on the pr-shape in the collider case. 

There is, however, no argument that enforces /ir and fip to be exactly equal, they are 
just expected to be of the same order of magnitude, given by the prompt photon's pp. 
In fact, /iR 7^ fiF automatically happens if one uses 'optimized' scales as mentioned 
in the introduction. A smaller renormalization scale /ir along with a larger factorization 
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scale jj,F can be expected to create a steeper slope of the theoretical result, since lowering 
IMr mainly increases the strong coupling constant a s , whereas the main effect of a larger 
fiF is to deplete the gluon distribution at larger x and to increase it at smaller x. In 
fact, the curves in Fig. 4 show that these effects are quite significant. The choices of, e.g., 
fiji = 0.3pT, 1±f = Pt or fi R = p T , fi F = 0.3pT do lead to about ±20% shape changes in 
the CDF region, respectively. Obviously, scale dependences are able to affect the px- slope 
of the NLO cross section, contradicting the corresponding conclusions drawn in fl8], [H[ 
which were derived assuming /i^ = /ip. All in all, scale changes seem to have a rather 
strong influence on the theoretical cross section even beyond the LO. This is in line with 



the observations in P7| , where it was shown that at small xt the scale dependence is 
only slightly reduced when going from LO to NLO, which renders it difficult to estimate 
the most appropriate scales. The scale dependence of the NLO cross section for prompt 
photon production indicates the importance of corrections of even higher order and sets 
severe limits on the accuracy of gluon determinations from these data. 



4 Combined Analysis of DIS and Prompt Photon 
Data 

In this section, we examine the question of whether the agreement between NLO cal- 
culation and the isolated prompt photon data can be further improved by adapting the 
proton's parton content, in particular its gluon density. For this purpose, we perform 
combined NLO fits to direct-7 production and DIS structure function data, the latter 
pinning down the quark densities, with different choices for the renormalization and fac- 
torization scales. In this way, we also investigate the uncertainty of the resulting gluon 
distribution originating from the scale dependence of the prompt photon cross section. 

Technically we proceed as follows: at the reference scale of Ql = 4 GeV 2 the gluon 
input is parametrized as 

xg{x,Q 2 ) = A g x a °{l -xf B (l+j g ^/x + 5 g x) . (3) 

This functional form is also used in the latest MRS analysis 0. For each given set of 
a g , j3g, 7 9 , and 5 g , a fit of the quark densities to F 2 P data is performed, where A g is fixed 
by the energy-momentum sum rule. Here the non-singlet quark densities u v = u — u, 
d v = d — d and A = d — u are, as in Q , for = 200 MeV directly adopted from the 
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MRS (A) global fit 0. For other values of A we employ the modified non-singlet sets of 
HH . This procedure guarantees a sufficiently accurate description of all observables testing 
mainly the flavor decomposition of the quark content and facilitates the fitting procedure 
since it reduces the number of free parameters. The sea quark input is parametrized as 

@ 

x(u + d)(x,Q 2 ) = A(.x a t(l-xft (l + 7 €v /^ + %r) , (4) 
xs(x,Ql) = -x(u + d)(x,Ql) . 

This input is evolved in the factorization scheme of @] (see section 3), calculating the 
charm contribution to F 2 via the LO Bethe-Heitler process using the LO gluon distribution 
and Alq of Q]. The five parameters in eq. (4) are fitted to the available F£ data of 
BCDMS 0, 0, NMC H, ZEUS 0, and Hi @ in the region where the structure 
function is sensitive to the sea quark and gluon densities, x < 0.3, and where higher twist 
contributions are expected to be small, Q 2 > 5 GeV 2 . As in [JJ] the BCDMS data are 
normalized down by 2%. At the present level of experimental accuracy, the HERA 



normalization uncertainties pll 22] can be disregarded for our purpose. Statistical and 



systematic errors are added quadratically for all data sets. 

We then use the complete set of parton distributions obtained from this _F 2 P fit for a 
fixed gluon shape to determine the Xdir.-r f° r the prompt photon data of WA70, R806/7, 
UA2 and CDF |TJ|, [15], already discussed in section 3, supplemented by the pp 



inclusive results of NA24 f§ and UA6 fT|], both at y/s = 2A GeV. Due to the very 



small point-to-point errors of the CDF data [18] their overall normalization uncertainty 
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is rather important. Therefore it is allowed to float in the fit with a contribution to Xd 
according to its experimental uncertainty of 10%, see ||. All other data sets, including 
P, [H| where separate normalization uncertainties are provided, are fixed at their nominal 
normalization. The whole procedure is repeated until finally Xtot = Xdis + Xdir.-y * s 
minimized. We have chosen this two-step approach of two 'nested' fits since it affords the 
least number of evaluations of x\i r . 7 which dominate the consumed computer time due to 
the complexity of the NLO prompt photon cross section calculation. In order to examine 
the uncertainty coming from the choice of the QCD scale parameter A, we perform fits 
for A^ = 200 MeV and 300 MeV. These values are representative of the present range 



of A found in analyses of DIS and related data 0, ^3], [44], [|5[]. In view of our findings 
for the scale dependence of the prompt photon cross section obtained in the last section, 
we repeat the fits for various combinations of the renormalization/f actor ization scales, 
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systematically scanning the range 0.3 pr < //r, Hf < 2.0 pr, see below for a detailed 
discussion. 



Fig. 5 displays the results of three of our fits, again compared to the data of pi [L4L 15 



18| . For clarity, the additional large-^T data sets of @, |Tj| are not shown in the figure. 



The upper part presents the fit using the 'standard' scales and of section 3, the 
lower plot depicts further examples of best possible fits, one for each A- value employed. 
The latter two fits happen to lead to the same CDF normalization factor of +10%. The 
input parameters and x 2_resu hs for these three representative fits are given in Table 1. 
It is obvious that we obtain good fits, with x 2 per data point slightly below 1 also in the 



prompt photon subset of the fitted data. The ISR and fixed-target data [p], Rl O, ITJ; 



are described about as well as in the reference calculation in section 3. The \ 2 for the 



CDF results JT8J seems to remain a bit high. We have not been able to reach a value 
below 25 for the 16 data points, however with more than 10 units contributed by just the 
measurement at = 48.9 GeV. The figure clearly demonstrates that no shape problem 



remains here. With respect to the UA2 results []15" |, the situation does only slightly 



improve; with present errors, however, the UA2 data do not impose strong constraints 
on the fits. On the other hand, it should be noted that a data set with about the same 
central values but much smaller errors would, in combination with the CDF results, very 
seriously challenge the NLO framework used in this paper, since we have been unable to 
improve the description sizeably even by artificially reducing the UA2 errors in special 
runs of our fits. For the larger value of A^-, the range of renormalization and factorization 
scales giving very good fits to the data is somewhat shifted towards larger values. As can 
be seen from Fig. 5 and the table, we do not find a significant difference in the quality of 
the fits for the two A-values. Taking into account the full scale dependence of the NLO 
cross section reduces the sensitivity of the prompt photon data to as as compared to, 
e.g., H [52] where 'optimized' scales were used. 



Our results demonstrate that presently published data on pp, pp — > 7X can be de- 
scribed quantitatively by NLO perturbative QCD. This finding is at variance with the 



results of a recent, partly comparable analysis of the CTEQ group f31| . In that paper, it is 
concluded that (at least) the j>r-shapes of the ISR and collider data cannot be satisfacto- 
rily fitted, unless an 'intrinsic' fc^-smearing is introduced, which is as large as < > ~ 4 
GeV for SppS and Tevatron conditions. It should be clear that we do not claim positive 
evidence for the absence of such a (somewhat counterintuitive) smearing, we just state 
that the data considered here can be accounted for without this assumption. It should 
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200 


200 
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l 1 Tt / Tim 


n 5 




n 7 


H>F/ FT 




n 7 


n 7 


a g 


-0.011 


-0.125 


-0.080 




5.715 


5.314 


5.024 


19 


-4.174 


-4.200 


-3.963 




5.217 


5.370 


4.931 


A 


U. ( oo 


U. 1 ID 


U. (o ( 




-0.156 


-0.148 


-0.135 


fk 


7.480 


7.641 


7.173 




0.642 


0.740 


0.843 


h 


1.521 


1.622 


1.080 


Xdis 


253.6 


257.3 


254.4 


Adir. 7 


59.1 


52.8 


56.5 


"X.tot 


312.7 


310.1 


310.9 



Table 1: Input parameters for eqs. (3) and (4) for three representative combined fits to 
direct-7 and F 2 P data, together with the resulting x 2 -values for the 60 prompt photon and 
294 DIS data points. 



be noted in this context, however, that a <kx> of about 1 GeV for fixed-target energies 
would spoil |46| the successful description of the WA70 data JTJ|. The consequences of our 
optimal 'pure' NLO fits to all direct photon data will be elucidated below. 



Let us note before some obvious differences between our study and ^TJ. With respect 
to data, we use the latest results of R806/7 fL4"| , whereas |31| includes previous superseded 
steeper results of this collaboration [jllj] . Taking into account the results of ||14j| , there 
is obviously no slope offset between NLO theory and ISR data, neither for 'standard' 
distributions, see Fig. 1 and the last reference of @, nor for our fits discussed above. In 
the fixed-target region, we do not include the recent E706 data [57] from pBe collisions 
in order to avoid any bias from a possible EMC-like effect on the nuclear gluon density. 
With respect to the theoretical treatment, in contrast to |H] we fully incorporate NLO 
fragmentation and the corresponding treatment of isolation into our analysis. The very 
small errors of the CDF data |l8j make a complete and consistent NLO analysis mandatory 
in order to arrive at any solid statements on the viability of the NLO description of prompt 
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photon production.^ 



Hence in the present situation it seems reasonable to us to maintain the 'pure' NLO 
perturbative QCD framework and to investigate the consequences of optimal fits to the 
direct-7 data in this scenario. In Fig. 6 the gluon distributions obtained in the three fits 
displayed above are compared at Q 2 = 20 GeV 2 to the corresponding result of the recent 
MRS (A') global study || where only fixed-target prompt-7 data were included in the 
analysis. Although data at different pr probe the parton distributions at different scales, 



the main effect of including the CDF results [[D| is readily seen in the figure. The gluon 
distribution is increased at the lower end of the x-range accessed by these measurements, 
x ~ 0.01, and is decreased around x ~ 0.15, thereby in between providing the steeper 
rise suggested by the data. At very small x, x ~ 10~ 3 , our gluon densities are larger 
than the one of MRS (A') and more similar in magnitude to the MRS(G) @ and GRV |§ 
parametrizations. 

Effects of such an enhanced gluon density at x ~ 10~ 2 should also show up in other 
processes which are sensitive to g(x, Q 2 ) in that x-region such as, e.g, bottom production 
in high-energy pp collisions at the Tevatron. Fig. 7 illustrates this fact. The cross section 
for single-inclusive 6-quark production, a(pp —> bX), expected to be mainly driven by 
gluon-gluon fusion, has been measured as a function of the minimal transverse momentum 
pmin Q £ 6-quark by the CDF and DO collaborations recently [HS , f|9| . The figure shows 



these data together with corresponding NLO calculations based on the work of |50], [Hj] 
for two of our fitted parton distribution sets and the MRS (A') parametrization For 
this illustration, we have chosen \i = \\Jpt + rn 2 for the renormalization and factorization 
scales and mj = 4.5 GeV for the bottom mass. Note that the theoretical uncertainties 
from variations of these parameters are rather substantial [jiS , [49| . 



Let us finally come back to the scale dependence of the fitted gluon densities. As 
already mentioned above, we have scanned the scale range 0.3 pr < fiR, Hf < 2.0 pt in our 
fits. The lower limit has been introduced in order to avoid low scales below about 1 GeV. 
We find very good combined fits to the complete set of prompt photon production and DIS 
structure function data given above for 0.3 (0.5) ~ fj,R, /ip /px ~ 0.7(1.0) for = 200 
MeV (300 MeV), respectively, and /ir ~ hf- The spread of the gluon distributions of all 

3 The input form used in the fits of |HJ has not been given in the paper. In the likely case that the 
CTEQ3 m input ansatz has also been employed in j3l[], another difference with respect to our study 
would be added, since we find that the inclusion of the -y^E-term i n (3) ; not present in ||, is important. 

4 We thank M. Stratmann for performing these calculations. 
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these best fits can serve as a measure for the theoretical uncertainty on xg(x, Q 2 ) induced 
by scale variations. Fig. 8 displays the resulting 'error band', normalized to the MRS (A') 
gluon distribution at Q 2 = 20 GeV 2 . Also shown are the uncertainties at four selected 
x-values arising from the experimental errors. They have been determined from fits for 
one scale combination, which are one unit in the total \ 2 higher than the corresponding 
best fit. In view of the controversial theoretical interpretation |30], |31| of the isolated 



high-energy collider data |15|, [X8f| , we have repeated the fits for the theoretical error band 
retaining only the 30 (non-isolated) fixed-target and ISR data points of |], H, [TTJ, [16| in 
the direct-7 part of the analyses. The resulting band is also presented in Fig. 8. 

It becomes obvious from the figure that the same pattern of deviations from the 
MRS(A') fit P|, namely an increase (decrease) at x around 0.01 (0.15), is present for all 
scale choices which allow for a good fit of the complete set of prompt photon data. Thus 
this effect is significant and cannot be removed by a 'better' choice of scales. Not too 
surprisingly, the fits restricted to the non-isolated data lead to gluon densities consid- 
erably closer to standard distributions. They are, taking into account the experimental 
uncertainties, in fact in good agreement with, e.g., the GRV gluon parametrization 0] 
employed as a standard choice in section 3. At x < 0.1, the experimental uncertainties 
and the theoretical error bands of the fits to all data are of about the same size. However, 
in the classical fixed-target regime, x > 0.3, mainly considered in || and many subse- 
quent analyses via 'optimized' scales f33|, the uncertainty of g(x,Q 2 ) is dominated by 



scale variations. For instance, at x ~ 0.4, the spread of the fitted gluon densities is as 
large as almost a factor of two. In this large-x region, shown separately in Fig. 9 also for 
a high scale typical for the production of particles with masses of a few hundred GeV, a 
comparison of recent parametrizations, e.g. of MRS(A',G) and GRV Q as displayed in 
Fig. 8, does not indicate the real uncertainty of the gluon distribution. 



5 Conclusions 



We have performed a detailed NLO perturbative QCD study of prompt photon production 
in pp and pp collisions with respect to its sensitivity to the proton's gluon content. In fact, 
our analysis is the first one that includes both the full NLO treatment of the fragmentation 
contribution and the complete NLO treatment of isolation, and the consideration of the 
full information available now from experiment. We have carefully studied the theoretical 



14 



uncertainties, especially those arising from the so far experimentally unknown parton-to- 
photon fragmentation functions and from scale variations. 

The dependence of the theoretical cross sections on the renormalization and factor- 
ization scales turns out to be the dominant source of uncertainty. Particularly at large-x, 
x > 0.3, the scale dependence very severely limits the accuracy of NLO gluon determina- 
tions from (fixed-target) prompt-7 production data. Present-day gluon parametrizations 
of CTEQ3, GRV, and MRS(A',G) are very similar in this region and their difference does 
not represent the real uncertainty of g(x,Q 2 ) here. The isolated direct-7 cross sections 
are also sizeably affected by the choice of these scales in magnitude and p^-shape. 

We have carried out combined global analyses of fixed-target and collider prompt- 
photon cross sections and DIS structure function data within the perturbative NLO frame- 
work. We find that such combined fits work very well and give a good description of all 
prompt photon data for various combinations of renormalization and factorization scales 
and values of the QCD scale parameter A. The gluon distributions resulting from such 
fits turn out to be larger at x 0.01 and smaller at x around 0.15 than those obtained 
in recent global analyses of parton distributions not taking into account the full set of 
prompt photon data. 
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Figure Captions 



Fig. 1: The 'default quantity' (a exp — o- th )/a t h vs. xj- = 2p T /y/s for the data of [1], 



14] , |15| , p~8|| as compared to the NLO theoretical cross section o~th, using the GRV 

6[ . The curves present 



parton distributions and photon fragmentation functions [ 
the shifts (cr th i — o~ t h) / 1 o~th, where a t h' denotes the theoretical cross section if the 
fragmentation contribution is neglected or if the fragmentation functions of [25] are 
used. 
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Fig. 2: Comparison of the fragmentation functions zD^(z, Q 2 )/a and zDJ(z, Q 2 )/a from 
ACFGP [H and GRV J26| at Q 2 = 10 and 100 GeV 2 . 



Fig. 3: NLO MS fragmentation piece for the ACFGP and GRV |2§ photon fragmen- 
tation functions at ISR and Tevatron energies, using the CDF isolation criterion in 
the latter case. Also shown are the contributions stemming from gluon-to-photon 
fragmentation only. All cross sections are normalized to the NLO direct piece. The 
parton distributions of || are used and the scales are Hr = Hf = 0.5 py. 

Fig. 4: Same as Fig. 1, but the lines displaying the shifts in the theoretical results if the 
renormalization and factorization scales are varied as indicated in the figure. 



Fig. 5: The data of [1], [14|, [T5|, |18] as compared to three NLO fits for the different choices 



of Ajyfg and of the scales indicated in the figure. The CDF data |1| are shown with 
their fitted normalizations. 



Fig. 6: Comparison of the gluon distributions obtained in the three fits displayed in Fig. 

MS 



5 with the MRS (A') [| gluon density. A denotes A^ in MeV. 



Fig. 7: The cross section for single-inclusive bottom production, a(pp — > bX), as a func- 
tion of the minimal transverse momentum p™ tn of the 6-quark as calculated in NLO 
perturbative QCD using the parton distributions indicated in the figure. Also shown 
are the (inclusive lepton) data of JJR |49[ . 



Fig. 8: The error bands of the fitted gluon densities due to scale variations for the fit to 
all data and for the fit retaining the (non-isolated) fixed-target |TJ, [5], [T(| and ISR 
[l4|j prompt photon data only, plotted via the deviation from the MRS (A') fit. Also 



displayed are the la uncertainties resulting from the experimental errors and two 
further 'standard' gluon parametrizations 0. [J] . 

Fig. 9: Same as Fig. 8, but for the absolute gluon distributions at large x at a typical 
scale probed by the direct-7 data (left) and at a very large scale relevant for the 
production of very heavy particles (right). Of the rather similar 'standard' gluons 
only the one of MRS (A') M is shown for comparison. 
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